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Abstract 

We studied the dynamics of isotropic-nematic transitions in liquid crystalline polymers by in- 
tegrating time-dependent Ginzburg-Landau equations. In a concentrated solution of rodlike poly- 
mers, the rotational diffusion constant D r of the polymer is severely suppressed by the geometrical 
constraints of the surrounding polymers, so that the rodlike molecules diffuse only along their 
rod directions. In the early stage of phase transition, the rodlike polymers with nearly parallel 
orientations assemble to form a nematic polydomain. This polydomain pattern with characteristic 
length £, grows with self-similarity in three dimensions (3D) over time with a £ ~ t 1 / 4 scaling law. 
In the late stage, the rotational diffusion becomes significant, leading a crossover of the growth 
exponent from 1/4 to 1/2. This crossover time is estimated to be of the order t ~ D~ x . We 
also examined time evolution of a pair of disclinations placed in a confined system, by solving the 
same time-dependent Ginzburg-Landau equations in two dimensions (2D). If the initial distance 
between the disclinations is shorter than some critical length, they approach and annihilate each 
other; however, at larger initial separations they are stabilized. 

PACS numbers: 64.70.mf, 61.30.Vx, 61.30.Dk, 61.30.Jf 



I. INTRODUCTION 



Liquid crystalline polymers (LCPs) are widely used in technology as high-performance 
fibers since they have high strength and are light in weight. Their mechanical strength 
increases when the high molecular weight polymers are orientated in the same direction 
[1-7J. Compared to low molecular weight liquid crystals (LMWLCs), LCP molecules do 
not align spontaneously. In this paper, we study the polydomain formation after quenching 
an LCP system from the isotropic to the nematic state. The phase transition dynamics 
in LCPs have been mostly studied on the basis of the conventional nematohydrodynamic 
equations [5HTU]. which were developed to describe the dynamics of LMWLCs [TT]. We are 
thus interested in the possible differences in the pattern evolution in LCPs and LMWLCs. 

In a concentrated solution, rodlike polymers entangle with one another, and hence, the 
surrounding polymers strongly suppress the rotational and perpendicular diffusions [3j. In 
the high concentration limit, each rodlike polymer moves only along its molecular axis and 
this parallel diffusion dominates the dynamics of the isotropic-nematic transition. 

If the rotational motion of the director is absent, the orientational order parameter be- 
haves as a conserved variable [Zj. It is well known that the mechanism of domain growth 
and the resultant growth exponent depend on whether its order parameter is preserved. For 
a system described by a single non-conserved order parameters such as magnetization, the 
domain pattern with a characteristic length £ grows in time as £ oc t 1//2 [T2] . On the contrary, 
when the order parameter is preserved, the domain growth obeys £ oc t 1 ^ 3 as observed in the 
phase separation of binary mixtures [T2] • The former is termed as "model A" and the latter 
as "model B" p]. 

In a typical LMWLC, each molecule can freely rotate to align parallel to the surround- 
ing molecules. The phase transition dynamics are well described by the time-dependent 
Ginzburg-Landau equation with a non-conserved tensorial order parameter [11J. In the late 
stage, the characteristic length of the polydomain pattern £ grows in time as £ oc t 1 ' 2 [TIHTB] . 
For LCPs, Shimada et al. studied the early stage of the phase transition and predicted a 
spinodal decomposition of the nematic order parameter [Tj. However, since their kinetic 
equations are linearized, the domain growth in the late stage could not be treated. In this 
study, we reformulate the free-energy functional and kinetic equations for rodlike polymers 
to include non-linear terms. This enables us to analyze the late stage behaviors of the phase 
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ordering, such as domain growth and defect motions. Similar kinetic models for mixtures of 
isotropic liquids and semi-flexible polymers have been proposed by several authors [17, 18J. 
However, these studies have focused on the phase separation in the mixtures and isotropic- 
nematic transitions have not yet been studied. The main aim of this paper is to elucidate the 
isotropic-nematic transition dynamics in concentrated solutions of LCPs and to determine 
the effects of small but finite polymer rotational diffusions on the system. 

This article is organized as follows: In section 2 we formulate free-energy functional and 
kinetic equations for the rodlike polymers in accordance with the method of Shimada et 
al. In section 3, we show the results of the numerical simulations on the isotropic-nematic 
transitions and discuss them. In section 4, we summarize our work. 

II. FREE-ENERGY FUNCTIONAL AND KINETIC EQUATIONS 

Similar to Shimada et al. [SI [7] , we derive the free-energy functional and kinetic equa- 
tions for a solution of LCPs. We reformulate their linearized equations to time- dependent 
Ginzburg-Landau equations [13] in order to study the late stage of the phase transition. 

Rodlike polymers of length d and width w (id < d) are considered. For a solution of the 
rodlike polymers, we introduce the free-energy functional J 7 for /(r, it), which represents 
the probability distribution of rods at position r where u is the unit vector along the rod 
direction. It consists of two parts as shown in Eq. (TO) 



J"o is the free-energy functional for ideal non-interacting polymers and is expressed as 



where T is the temperature, fee is the Boltzmann constant, and vo(= w 2 d) is the volume of 
the rodlike polymer. The second term in Eq. is the interaction part of the free-energy 
functional and is expressed as 



J 7 — -7"o + -F int- 



(1) 






xf(r,u)f(r',u). 



(3) 
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W(r, u, r', u') represents the excluded volume interaction between two rodlike polymers 
(r, it) and (r', it'), and is defined by 

{a (when two polymers intersect) 
(4) 
(otherwise), 

where a is the interaction parameter that has the dimension of volume and is estimated to 
be wd 2 [51 119]. We neglect the interaction between the polymers and solvent; therefore, the 
isotropic-nematic transition originates purely from the configurational entropy of the rods 
[TT| [20]. In other words, the solution is an athermal system, in which temperature changes 
play no role in the phase behaviors. Above a critical concentration, the solution exhibits a 
liquid crystalline phase. This free-energy functional is applicable not only to LCP solutions 
but also to suspensions of rigid rods, such as the tobacco mosaic virus (TMV) [2TH23] and 
carbon nanotubes |24j. 

We define two order parameters as 

0(r) = v J duf(r,u), (5) 

Qij(r) = v j duf(r,u) (^UiUj - \ , (6) 

where J du represents a solid angle integration. <p is the concentration of the polymers, and 
Qij is the orientational order per volume. It is noted that Qij vanishes with vanishing 0, even 
when the polymers are orient at ionally ordered. Although it is more natural to use Qij/<f) as 
an order parameter, we evaluate the free-energy functional with <fr and Qij for simplicity. 
The distribution function can be expanded for and Qij as shown in Eq. ^ 

f(r, u ) = ^~ o j^) + y% (r) [uiu 3 - \8iA } • (7) 

Hereafter, the repeated suffixes i and j indicate summation over i, j = x, y, z. We substitute 
Eq. ([7]) into Eq. ([2]), and integrate it only over u using the isotropic approximation (see 
Appendix A). After some calculations, the free energy for ideal polymers J-'o is expanded to 
include the fourth order of Q^ as shown in Eq. (pi) 



c 2 d \ 

^^^(QijQji) + -^QijQjkQklQlij , (8) 



where the numerical constants are ao = 15/2, b = 225/14, c = 375/28, and d = 375/7. 

Using the isotropic approximation, we obtain the interaction part of the free energy [Eq. 
(pi)] in reciprocal qr-space as follows: 
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+^^d 2 {q k q k Q ij (q)Q ij (-q) + 4|g i Q ij (<?)| 2 } 



(9) 



where X(q) is the Fourier qr-component of the variable X(r) in reciprocal space. The expres- 
sion of W in reciprocal space is denoted in Appendix B. With reverse Fourier transformation, 
we finally obtain T^t in real space as 

k B T 

J~ int. — 




+ y|V0| 2 + y|V 2 0| 2 + KoVrfVjQij 

+ y \^iQjk\ 2 + I Y |ViO«| 2 J, (10) 

where x = 7ra/(2uo), ai = 157ra/(32t>o), Co = — 7ra<i 2 /(72t>o), Ci = ll7ra<i 4 /(57600t'o), Kq = 
77iad 2 /(768v ), K x = 15vrad 2 /(2688w ), K 2 = 15tt ad 2 /(672v ), and V* represents 9/9r< 
(z = x,y,z). In our model, the system is in the isotropic state for < 0i(= 976t> /637ra) 
and in the nematic state for <fi > (f>2(— 16fo/vra). When <fii < <fi < 2 , both the phases 
coexist. 

It is noted that the gradient terms of are expanded up to the fourth order of q because 
the coefficient of q 2 \4>\ 2 is negative in Eq. This negative coefficient implies that the 

density modulations have a periodicity of 2n(— Ci/Co) 1 / 2 , and its contribution on the phase 
ordering is small as discussed in Appendix C. However, actually, our numerical simulations 
do not show such a modulated pattern in our concentration range. 

Next, we introduce the auxiliary fields /x and Hij as 

jj- = v o + (uiUj - i<y J I . (11) 

With isotropic approximation, the variation of the free-energy functional is given by 

5f= J driSfa + SQijHij). (12) 



From this equation, the following expressions are derived, 

5J 7 „ 5J 7 



H, 



'.i 



5Q l 



(13) 



in which \x and can be interpreted as the chemical potential of <fi and the molecular force 
field of Qij, respectively. We should note that our definition of is different from the 
conventional definition in terms of its sign. 

Using the free-energy functional J 7 , the Fokker-Planck equation (see Eq. (2.1) of Ref. 
[7]) is rewritten as 



0_ 
di 



f 



knT 



V • {D\ 



uu 



„ -, 5 J- 5 J- 

D ± (I-uu)}V— + D T K 2 



(14) 



5f 5/. 

Here, <fi is the average concentration of the rodlike molecules and 1Z = u x (d/du) is the 
rotation operator [3]. Dm and D± are diffusion constants for parallel and perpendicular 
motions to the rod direction u, respectively, and D r is the coefficient for the rotational 
motion. For a dilute solution, the Kirkwood theory estimates the diffusion constants as Djj = 
k B Tln(d/w) /{2irr] s d), D° ± = k B Tln(d/w)/(47cr] s d), and D r ° = 3k B T{ln(d/w) - 7} / ^(vr^rf 3 ) . 
Here, rj s is the solvent viscosity and 7 is Euler's constant [3j. As already noted, in the 
high concentration limit, each rodlike polymer moves only along its molecular axis, namely, 
Dm ^> D r d 2 ,D±. Hereafter, we use the kinetic coefficients Lx = Dx4>vo/k B T , where X 
represents || , _L, and r. 



From Eq. (14), we express the kinetic equations for and Q~ with isotropic approxima- 



tions as follows 
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(16) 



The off-diagonal coefficients in Eqs.(15) and (16) satisfy the Onsager's reciprocal relation- 



ship, and in Eq. (15), we omitted V 2 Hij5ij/3. The time derivative of the free-energy 
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functional is given by 

d 
dt 



- / dr 



dt r 
(V./i) 2 
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(17) 



In the second line of Eq. (17), we ignore the influence of and Qij from outside the system. 
The conditions Lii > Lj_ > and L r > guarantee that the integrand of Eq. (17) is positive, 
resulting in dF /dt < 0. 

We normalize space and time by d and to — d 2 / (-^110)- I n an aqueous suspension of TMV, 
we estimate that d = 320 nm and w = 20 nm. Assuming T = 300 K and r] s = 0.9 mPa • s, 
the diffusion constant is D\\ = 6.3 /im 2 /s, and we thus obtain to — 1O2.40 -1 ms. We integrate 
the coupled equations in the lattice space with the explicit Euler method. In order to save 
computational costs, the spatial Ax and temporal At increments are varied according to the 
state point. In all the following simulations, we employ periodic boundary conditions. 



III. RESULTS AND DISCUSSION 
A. Spinodal decomposition 

First, we study the isotropic- nematic transition in a concentrated solution of rodlike 
polymers. We perform 3D simulations and set a = 16t>o, which corresponds to d/w = 16. 
The spatial and temporal increments are Ax = OAd and At = 0.016to, respectively. As an 
initial condition, we set = 0.5, which is larger than the critical concentration 2 - The 
spatial averages of all the components of are set to zero and we add random noise to 



them. Here, we set L± = L r d 2 = in Eqs. (15) and (16), so that the phase transition 
proceeds in the rodlike polymers only via the diffusions along their axes. 

Shimada et al. reported spinodal decomposition-like growth of the nematic order param- 
eter in the early stage of phase ordering [7] . They also claimed that the domain growth can 
be separated into three modes, i.e.,, splay, twist, and bend. Accordingly, we decompose the 
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structure factor of Qij into 



&pi(g) 



Stws(q) — 2\a-iQijbj\ 2 + 2\diQijC 



'3 I ' 



Sbnd(<?) = \biQijbj - CiQijCjl 2 + AlbiQijCj 



(18) 
(19) 
(20) 



Here, a(= q/\q\) is the unit vector toward the wave vector q, and b and c are also unit 
vectors, which are orthogonal to a and each other. 

Figure [TJ shows the decomposed structure factors of Qij at t = 2.56to and 5.12t - Their 
shapes are similar to those found in the spinodal decomposition of phase separations. 
Namely, each structure factor has a peak at an intermediate wave number and vanishes 
for q ~ 0. Neglecting the higher order terms, we obtain the early-stage growth rates of 



the three modes in Eqs. (18)-(20), as denoted in Appendix C. The positions of the peaks, 



predicted by the linearized analyses, are marked by arrows in Fig. [TJ The simulation results 
are consistent with the linearized theory and the splay mode develops more slowly than the 
other two modes. 
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FIG. 1: (color online) Decomposed structure factors of the nematic order parameter in isotropic- 
nematic phase transitions at t = 2.56£o (empty symbols) and 5.12io (full symbols). Black circles, 
red squares, and blue triangles represent the splay, Twist, and bend modes, respectively. Peak 
positions predicted by the linearized analysis are marked by the arrows. 

The linearized analysis indicates that the fluctuation of <fi and the splay mode are coupled 
to each other. However, our simulated structure factor does not show an appreciable 
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peak (data not shown). We consider that this is an artifact of our numerical simulation, 
because the complex spatial operators in the kinetic equations are difficult to deal with 
precisely. We need to improve the numerical scheme to study the spatial distribution of </> 
more quantitatively. 

B. Growth of polydomain for L v = 

Next, we study the temporal evolution of the polydomain pattern in 3D. We set a = 16v o, 
Ax = d, and At = 0.05to- The initial condition is = 0.5, which is larger than 02, and 
Qij = with the random noise. Figure [2^ a) shows the temporal evolution of the pattern 
of Q 2 xy in an xy-plane (z = 0). The corresponding director field at t = 8000to in the same 
xy-p\ane is shown in Fig. ^b). A polydomain pattern is formed and it coarsens in time. 
Figure [3] schematically explains how the isotropic-nematic transition takes place without 
molecular rotation. If the rotational motion is allowed, the rodlike polymers rotate to align 
with fixed positions as shown in Fig. [3]^a). On the other hand, when the rotational motion 
is severely suppressed, the rodlike polymers that are initially nearly parallel to each other 
assemble to form a small grain via diffusion along their axes. The assembled grains form 
a mosaic pattern and many defects remain as shown in Fig. ^b). In the latter case, the 
nematic order parameter is conserved in the whole system. 

The defects in Fig. [2] are at the intersections of disclination lines in the xy-plane (z — 0), 
entangled in three dimensions. In 2D Schlieren textures, the number of bright brushes 
forming a defect core is given by 4|m|, where m is the topological strength of the defect. 
We observed that most defects have two brushes in Fig. |2](a). As in other nematic states of 
LMWLCs, disclination lines of m = ±1/2 are formed more frequently than other types of 
defects. 

Since the molecular rotation is severely suppressed in LCPs, its coarsening mechanism is 
very different from that in LMWLCs. Even after the early stage, the scalar nematic order 
parameter, which is given by Q° = (2QijQji/3y^ 2 , remains inhomogeneous; usually, Q° is 
smaller than the equilibrium nematic order Q eq near the defects. Since the inhomogeneity 
of Q° affects the structure factor in the high g-range, we calculate the structure factor of a 
normalized order parameter, = Qij/Q°, to determine the evolution of the polydomain 
pattern. This normalization corresponds to a binarization method for phase separation [25J . 
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FIG. 2: (a) Snapshots of Q 2 xy in an xy plane (z = 0) at t/t = 400, 1000, 2000, and 8000 in the 
isotropic-nematic transition without rotational diffusion, (b) The corresponding director field of 
the domain pattern in the same xy plane at t = 8000to- The lines represent the principle axes of 
the tensorial order parameters. 

In Fig. [i](a), we plot the temporal change of the total structure factor S tot (q) = \Qij{q)\ 2 , 
in which Qij{q) refers to the Fourier transform of Qij{r). The structure factor is not decom- 



posed into the three modes given by Eqs. (18)-(20). It is shown that the peak position of the 



structure factor shifts toward q = and the peak height develops with time. These features 
are similar to those in the late stage of phase separation [121 12S] . It is known that the struc- 
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(a) 




FIG. 3: (color online) Process of isotropic-nematic transition in a solution of rodlike polymers, 
(a) When the rotational motion is allowed, the polymers rotate to align with fixed positions as 
in LMWLCs. (b) In the absence of rotational motion, the rodlike molecules diffuse only along 
their rod direction to form small grains of the nematic phase. In the center column, red and blue 
segments represent initial (left) and ordered (right) configurations, respectively. Arrows of broken 
lines indicate the motions of the rodlike polymers. 

ture factor scaled by the characteristic wave number collapses into a master curve at isotropic 
phase separation. We replot the scaled structure factor, Stot(<?) = (9(^)) 3 Stot(<?/(<?(£))), 111 
Fig. [41(b) and the characteristic wave number is defined as 

/ dqS(q,t) tot \q\ 



<?(*)> 



(21) 



/ dqS tot (q,t) 

Figure |4^b) shows that the dynamic scaling law holds fairly well in the isotropic-nematic 
transition of rodlike polymers. 

In the high wave number range (q > (q(t))), the structure factor decays as S to t{q) ~ q 5 - 
In phase separation, a decay S(q) ~ q~ 4 , termed as Porod's law, which originates from the 
scattering of 2D interfaces in a 3D matrix, has been observed [T2] . 

The g~ 4 -tail was also reported in 2D simulations of a nematic phase [151 12SJ. The g~ 5 -tail 
observed in LCPs is considerably different from these tails. In similar systems, Bray studied 
phase transitions described by a conserved N- vector order parameter [27], and showed that 
the structure factor exhibited a q~( N+3 ^-tai\ in 3D. From his work, it follows that a g _5 -tail 
is obtained for a system of iV = 2 or XY model, in which many line defects are formed. We 
consider that the g~ 5 -tail observed in LCPs represents scattering from the entangled one 
dimensional (ID) disclination lines in the 3D matrix. 
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FIG. 4: (color online) (a) Temporal changes in the total structure factor of the nematic order 
parameter in isotropic-nematic transitions without rotational diffusion, (b) Structure factor Stot(<?) 
scaled by the characteristic wave number (q(t)) of the domain pattern. Stot(q) shows the domain 
growth with self similarity and the Porod law Stot ~ Q~ 5 is observed in the high q-regime. 

The self-similarity in the scaled structure factors indicates that the polydomain growth 
is characterized by only one characteristic length scale. We plot the temporal change of the 
characteristic polydomain size, which is defined as £(t) = 2n/ (q(t)) in Fig. [5j After the early 
stage, the characteristic size develops with time as £(t) ~ t a with a = 1/4. This exponent 
is smaller than those in phase separation (a = 1/3) and isotropic-nematic transition of 
LMWLCs (a = 1/2) [HHT6]. Interestingly, the exponent is the same as that of a system 
described by a conserved XY model [27| 128] . This coincidence and the same Porod's tail 
value imply profound similarities between LCPs and the conserved XY model. Therefore, 
the analysis of the conserved XY model might be helpful to understand the coarsening 
mechanism in LCPs. However, there is also an important difference between the two; while 
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most defects have m = ±1/2 topological strengths in the nematic state of rodlike polymers, 
the topological strengths in the XY model are m — ±1. Further studies are needed to 
clarify the similarities and differences between them. 




FIG. 5: (color online) Time evolutions of the characteristic domain length £{t). Black and red 
symbols represent those for L± = and L± = Ln, respectively. Rotational diffusion is set to 
L r = 0. In both the cases, characteristic lengths grow with time according to £(t) ~ t 1//4 . 



We have also studied the effect of the perpendicular diffusion on the polydomain growth. 



We set L± = Ln and L r = 0, where the off-diagonal terms of the kinetic equations (15) and 



(16) vanish. Since the rotational diffusion is not included, the tensorial order parameter 
Qij is still conserved. The numerical simulations show that the dynamic scaling law holds 
and the growth exponent is also given by a = 1/4 as shown in Fig. [5] (red squares). It is 
indicated that this growth exponent is not characteristic of the parallel diffusion, but stems 
from the nature of the preserved order parameter. The characteristic length for L± = Ln 
grows faster with time than for L± = by a factor of approximately 1.6. 



C. Growth of polydomain for L r > 

Solutions of LCPs have a very small but finite rotational diffusion coefficient D T . Thus, 
we expect that the polydomain growth will be affected by the rotational diffusion in the late 
stage of phase ordering. Figure [6] shows the temporal change of the scaled structure factor, 
Stot(<z)- Here, we set L± = and L T d 2 /L\\ = 0.0625 and the other parameters are the same 
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as those for L r = L± = 0. In the early stage, the structure factor has the same features as 
those in Fig. [4J namely, S tot (q) is very small at q ~ and it has a peak at an intermediate 
wave number. With time, the structure factor in the lower (/-range develops and in the late 
stage, S to t(q) has the Ornstein-Zernike form S(q) ~ (1 + g 2 £ 2 ) -1 , which is also observed 
in LMWLCs. The dynamic scaling law does not hold during the whole phase transition 
process and the growth of S to t{q) in the lower g-range is attributed to the rotational motion 
of rodlike polymers. 




FIG. 6: (color online) Time development of the scaled structure factor, <Stot(<7)> in the isotropic- 
nematic transition with rotational diffusion, L r d 2 = 0.0625L||. Gradual evolution in the low q- 
regime is observed after the crossover time t cr = lOOio- 

In Fig. [7| we show the time evolution of i for a number of L r 's. Although the physical 
meaning of £, especially in the crossover period (see below), is not clear, it is still a useful 
measure for the pattern growth. In the early stage, it does not change with time. This 
steady length corresponds to the spinodal decompositiondike growth of the nematic order 
parameter. After the early stage, the domain length evolves obeying £ ~ t 1 / 4 as in the case 
of L r = 0. Figure [7] indicates that the growth exponent changes from a — 1/4 to a — 1/2, 
which is the same as in the case of phase transition of LMWLCs, during a crossover period. 
Fig. [7] also indicates that the crossover depends on the rotational diffusion constant. As L r 
increases, the crossover is observed at earlier times. 

Each rodlike polymer moves along a tube surrounded by other tubes. The tubes are not 
necessarily straight and they disappear in a certain period of time, New tubes are continually 
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FIG. 7: (color online) Temporal changes in the characteristic domain lengths in the isotropic- 
nematic transitions in LCPs. Rotational diffusion constant is varied as L r d 2 /L\\ = 0, 0.0125, 0.0625, 
and 0.625. The occurrence of crossovers is indicated by arrows , and their positions are given by 
U = 6.25/A- 

created as the surrounding molecules fluctuate. As a result, the polymer gradually loses its 
original orientation [29]. We estimate that the crossover time t cr is of the order of the 
characteristic rotational time D~ l . After the crossover, the orientational order parameter 
is no longer conserved. In Fig. [7j arrows mark the corresponding crossover times and 
we assume t cr = 6.25t L\\/(L,.d 2 ). Although it is difficult to determine exactly when the 
exponent changes from a = l/4toa = l/2, the values marked by the arrows appear to be 
consistent with this interpretation. 

In Fig. |8j we show the time evolution of the Schlieren pattern Q 2 y in an x|/-plane for 
L v d 2 /L\\ = 0.625 and the crossover time is estimated as t CI = 10to- In an LCP solution 
with a finite L T , the director field rotates slowly, but freely, to adjust to the surrounding 
molecules. 

D. Defect motion 

We study the finite-size effects on the stability of defects, mediated by the elastic field 
of the nematic phase for a defect pair of anti-signed topological charges. As there is an 
attractive interaction between the charges, the defects approach and annihilate each other. 
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FIG. 8: Snapshots of Q 2 xy in the xy plane after the crossover. Rotational diffusion constant is 
L r d 2 = 0.625L.L. 

It was reported that the separation of R between the defects decreases to zero with time as 
R ~ {ta — t) 13 with (3 = 1/2, where t a is the annihilation time [2E1 EHl EH] • After annihilation, 
the director field relaxes to a homogeneous state in order to release the elastic energy in 
LMWLCs. However, this argument is not applicable to LCPs, since the kinetic mechanism 
is different. 

In the simulations, a = 16vq, = 0.5, Ax = 0.2d, and At = 0.0002to, and we set L± = 
and L T = to consider only parallel diffusions, i.e., those along the polymer molecular axis. 
We place a pair of anti-signed defects (m = ±1/2) in a small square box and carry out 2D 
simulations to mimic two parallel disclination lines with topological strengths of m = ±1/2. 

Initially, the spatial distribution of Qij is set to Qij = 3Q eq (uiUj — 5^/3) /2 with 



Here, 9+{r) is the angle between (r — r+) and (r_ — r + ), and vice versa. r + and r_ are 
the defect positions of m — 1/2 and m = —1/2, respectively [30]. At t — 0, (r + + r_)/2 
is at the center of the box and r + — r is along the x-axis. This configuration is not the 



u(r) = (cos{(^ + - 6„)/2},sm{(6 + - 0_)/2},O). 



(22) 
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equilibrium structure of the defect positions, because the spatial variations of the scalar 
nematic order parameter and the non-linear terms are neglected. The absence of these 
terms do not affect our results, since the structure relaxes quickly before defect motions are 
excited. In order to avoid numerical artifacts near the boundary wall, we employ the periodic 
boundary condition; this may be inappropriate to study the realistic confinement effects on 
the defect motion; however, as there are only two defects in the system, the simulations give 
valuable insights into the stability of the defects. 

Figure [9] shows time evolution of the director fields in Qij. In Fig. |9]^a), the initial defect 
separation is R(0) = \r + — r_| = 12 Ax and the box is H 2 , where H = 48Ax. The two defects 
approach and annihilate each other at t = 153. 3to- Contrary to the cases in LMWLCs, the 
director field remains distorted even after a long annealing time. This is because the y- 
component of the director field is preserved. In other words, the rodlike polymers, which 
are initially oriented along the y-axis, remain permanently aligned with the axis if L T = 0. 
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FIG. 9: Director fields around a pair of anti-signed defects at t/to = 8, 80, and 320. The initial 
separation between the defects is (a) R(0) = 12Ax and (b) R(0) = 16 Ax. The square box size 
is H = 48Ax and rotational diffusion is not allowed. Owing to the conservation of the order 
parameter, the director field remains deformed even after the defects are annihilated. At large 
separations, the defect structure is stabilized. 
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The free energy of the final distorted state in Fig. |9Ja) is lower than that of a uniform 
nematic state with the same average order parameter Q^. Here, it is important that 
differs from the equilibrium value Q^J. The final distorted state is determined by the balance 



between the local and non-local terms in the free-energy functional [Eqs. (J8J) and (10)]. For 
LMWLCs, the director field in the equilibrium state can optimize both parts of the free 
energy, such that Qij = Q e ^. On the other hand, as the lowering of the local part has to 
induce the deformation of the director field in LCPs, it is reasonable to assume that the 
elastic distortion remains even in the final state. 



Figure 10 plots the free energy at t — 800t calculated by Eqs. (|8j) and (10) versus the 
initial separation in a square box of H = 32 Ax. With increasing R(0), the free energy 
is increased for R(0) < 11 Ax. This increase enhances the elastic distortion without the 
formation of defects, as observed in Fig. [9j A kink in the free energy is observed around 



R = 11.5Ax. In Fig. 10, the free-energy difference from a reference {J r (llAx)+J : '(12Ax)}/2 
is plotted. This kink represents the critical separation R t , above which the defects do not 
annihilate each other and remain even after a long annealing time (t = 800£o)- The snapshots 
of the stabilized director field for R(0) = 16Ax in the box of H = 48Ax, are shown in Fig. 
(9^b). These stable defects also stem from the conservation of the order parameter. As the 
initial separation increases, the amount of the rodlike polymers oriented along the y-axis 
increases, so that the elastic field is appreciably distorted. Above the critical separation, 
the formation of local singular points (defects) are preferable to gradual distortion without 
defects. 

The critical separation depends on the system size. Figure pTj^a) shows the defect po- 
sitions as a function of time. We fix the initial separation to R = 12Ax and vary the 
system size by H/Ax = 32, 48, and 64. In the largest system, the defects approach faster 
and as H decreases, the defect motion becomes slower and the resultant annihilation time 
is retarded. For H = 32Ax, the defects initially experience a small shift at early times 
(t < 10to) and then hardly move. This dependence on the system size is unique to LCPs, 
and is not observed in LMWLCs. When the system size is large, there is a lot of room for 



the incompatible rodlike polymers to diffuse. In the inset of FigjlOj we show the dependence 
of the threshold Rt, on the average order parameter difference, SQ = f dr(Q yy — Q xx )/H 2 . 
In the initial configuration, the rodlike polymers along the y-axis are localized in between 
the two defects; therefore, the total amount of polymers is expected to be proportional to 
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FIG. 10: Plot of the free-energy difference at t = 800io with respect to the initial separation of the 
defect pair. The box size is fixed at H = 32Ax and the jump around R(0) — 11.5Ax represents 
the threshold Rt, above which defects remain even after a long annealing time. Inset: reduced 
threshold R t Ax/H 2 vs. 5Q. 

the defect separation, i.e.,, f dr{Q yy — Q xx ) oc R. When \SQ\ is smaller than a critical value, 
the solution cannot relax to a homogeneous nematic phase without defects and the simple 
scaling relation for the critical defect separation is R t /H 2 oc SQ. 

This size dependence is similar to that in systems described by a single scalar order 
parameter. Here, we consider a system whose free energy has two minima below its critical 
point. We assume that a droplet of one phase is placed in a matrix of the other phase and if 
the order parameter is not conserved, the droplet will be adsorbed in the matrix phase, as in 
magnetism. When the order parameter is conserved, as in phase separation, the droplet can 
stably exist in a confined system. In its steady state, the radius of the droplet is determined 
by the average volume fraction of the components, and the concentration of the matrix phase 
is slightly supersaturated compared to the equilibrium concentration. This supersaturation 
is related to the interface tension given by the Gibbs-Duhem relationship (12]. As the 
volume of the confined box increases with a fixed droplet radius, the droplet evaporates and 
the system becomes homogeneous. This is because supersaturation decreases with increasing 
box size and the resultant critical droplet size is increased. In LMWLCs, the supersaturation 
relaxes locally and quickly to the equilibrium state. On the other hand, in LCPs, the system 
is "supersaturated" from the equilibrium state and the large supersaturation leads to the 
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FIG. 11: (color online) (a) Time evolutions of the positions to the right (m = 1/2) and left 
(m = —1/2). Initial separation is R(0) = 12Ax and the box size is varied as H/ Ax = 32, 48, and 
64. When H = 32Ax, the defect positions are almost fixed after a small shift in the early stage, 
(b) Separation of the defect pair as a function of t a — t, where t a is the annihilation time. Initial 
separation and the box size are R(0) = 12Ax and H = 48Ax, respectively. The curve is fitted 
with the function (t a — t)^ , where i a = 153. 3to and (3 = 0.256. 



deformation of the director field. 

In Fig. |ll|b), we replot the defect separation as a function of reduced time, t a — t. The 
initial separation is R = 12Ax in a square box of H = 48Ax, and the defects annihilate 
each other at t = t a (= 153. 3t ). We estimate the exponent (3 = 0.256 in R ~ (t a — t)@ by 
fitting the curve. This is considerably smaller than = 0.5 for LMWLCs and close to the 
growth exponent a. Note that the annihilation exponent /3 is also of the same value as the 
growth exponent a in LMWLCs; however, the mechanism is not clearly understood and we 
have not concluded whether this exponent is universal for LCPs. Interestingly, Fig. [Tl^a) 
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also suggests that the defect motion before the annihilation becomes asymmetric; the defect 
for m = 1/2 moves faster than the defect for m = —1/2. Asymmetric motions of defects are 
also observed in LMWLCs with hydrodynamic interactions [50]; however, these interactions 
are absent in our model. 

IV. CONCLUSION 

We studied the isotropic-nematic transition in liquid crystalline polymers by integrating 
the time-dependent Ginzburg-Landau equations for the compositional order parameter 
and the orientational order parameter Qij. The kinetic coefficients are evaluated using the 
Fokker-Planck equation [7J . This approach ensures that the rodlike polymers diffuse only in 
the direction parallel to their molecular axis. 

Even if rotational motion is not allowed, the polymers being nearly parallel to each other 
assemble to form a nematic grain in the early stage of the phase transition. Since the director 
field is randomly oriented before the quenching of the isotopic phase, a polydomain structure 
is formed and many defects remain at the grain boundaries. The polydomain growth then 
exhibits self-similarity and the growth exponent is a = 1/4, which is considerably smaller 
than that for nematic liquid crystals of low molecular weight molecules. This small exponent 
is similar to those found in systems with conserved vector order parameters [27] ■ Here, the 
structure factor of the nematic order parameter has a peak at an intermediate wave number 
and nearly vanishes at q ~ 0. Hence, the orientational order parameter is preserved, contrary 
to that in LMWLCs. 

We have also shown that small but finite rotational diffusions can dominate the dynamics 
after a crossover time t cr . After the crossover, the growth exponent changes to a = 1/2, 
which is same as that for LMWLCs, and the polydomain pattern and the structure factor 
are also similar to those of LMWLCs. We estimated the crossover time as t CI = 1/D T , 
which enables rodlike polymers to rotate and the system behaves as a normal nematic liquid 
crystal, as in LMWLCs. This estimation qualitatively explains our numerical results. 

We have also shown that defect motion is strongly influenced by the conservation of the 
nematic order parameter. In LCPs, the defects can be stabilized in a confined system and 
the stability depends on the box size. When the stability is lowered by increasing the box 
size, a pair of anti-symmetric defects annihilate each other. The director field is distorted 
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after the annihilation of the defect pair and the defect annihilation obeys the power law, 
i.e., R~ {t^-ty with /3 = 0.25. 

There are many experimental studies on the phase transition of LCPs. However, most 
have been analyzed with the same kinetic equation used for LMWLCs. Before the crossover 
is reached, the molecule diffuses by a crossover length 4r = (D\\/D T )V 2 . For a semi-dilute 
solution of rodlike polymers, the crossover length £ a is comparable to the molecular length 
d. In order to examine the crossover in the phase transitions in LCPs, one must carefully 
probe the structure factor around q ~ 27r/d at very early times. In a concentrated solution, 
the rotational diffusion constant is approximated by D r = (3D®(<f)d 2 /w 2 )~ 2 , where /3 is the 
numerical factor [291 E2] . Hence, the crossover length is increased with the average volume 
fraction as £ CT ~ <pd 3 /w 2 . In a melt of LCPs or a highly concentrated solution of rodlike 
polymers, the crossover might be experimentally accessible, and we hope that our numerical 
study will stimulate detailed experimental observations in these systems in the near future. 
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Appendix A: Isotropic approximation 



When we derive the free-energy functional (Eq. (10)) and kinetic equations [Eqs.(15) and 



(16)], we employ the isotropic average approximation [7]. The isotropic average of A{u) is 
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defined by 



47T 

Then, the following formulas are obtained: 



(A) u = ^ I duA(u). (Al) 



(ui) u = (uiUjU k ) u = 0, (A2) 
(uiUj)u = ^ij, (A3) 
(uiUjU k ui) u = ^(5 ij 5 kl + 8 ik Sji + 5 u 5 jk ). (A4) 



Furthermore, we obtained 



(|u x u'\) u> = * (A5) 

7T 

(\u x u'^u'^w = - — (uiUj - 3%). (A6) 



In deriving Eq. (10), we need higher-order moments of u such as (|w x u'\u' i u' j u' k u' l } Ul which 



are given in Ref [7j. 

Appendix B: Excluded volume interaction between rodlike polymers 

In Eq. (|3]), W (r , u; r' , u') represents the excluded volume interaction potential between 
rodlike polymers (r,u) and (r',u f ). In Fourier fe-space, it is expressed by [HI 



W(k,u,u') = J drW{r,u,Q,u')e' 



kr 



sin(K ■ r) sin(K ■ u') _ , 

= 2a u x u'\ — ^ -— ^ Bl 

1 1 K u K u' y ' 

where a is the interaction parameter that is estimated as wd 2 . K is defined as kd/2, and d 
and w are the length and width of the rodlike molecules, respectively. We expanded it up 
to the order of k 4 using the following expansion: 



sm(K-u) (K-u) 2 (Ku 



K u 6 120 

This approximation is allowed for k < 1 jd. 



1 - r + 1on + • ■ ■ • (B2) 
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Appendix C: Linear analysis 



Here, we derive the growth rates of the order parameters in the early stage of the isotropic- 
nematic transition. We assume Lj_ = L r d 2 = in the high concentration limit. The 
deviations of the order parameters from the initial conditions are so small that the chemical 
potential and molecular field are linearized as 

fi = (^ + X + C q 2 + C x tA 50 + K^SQa, (CI) 

Hij = ( ^ - aij SQij + K qiqjS(j) + Kiq 2 5Qij + ^K 2 {qiq^Qkj + qjqk^Qki) , (C2) 

where 5(f) = <fi — (p and 5Qij = Qij. 

As noted in the main text, the orientational order parameter can be decomposed into the 
three modes splay, twist, and bend. According to Shimada et ai, the kinetic equations are 
also linearized as 

d ~ 

—5(f) = -Y^{q)8(j) - T^ s (q)5Q spX (C3) 
d ~ 

T^Qspi = -T s ^(q)8(f) - T StS (q)5Q sp i (C4) 
d ~ 

^5<5tws = -r t , t (g)<5Qtws, (C5) 
d ~ 

■7^SQ hnd = -r b ,b(<?)<5<2bnd, (C6) 

where 5Q sp i, SQ tws , and SQhnd are the decomposed splay, twist, and bend modes of the 
tensorial order parameter in Fourier space, respectively. The decomposed structure 
factors in Sec. Ill A are obtained from them. For example, 5Q sp \ stands for aiQ^aj in Eq. 



(18). The coefficients are given by 

1 

3 

f 
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r+Jq) = \h\<?{A + [C + ^K )q' + C iq A }, (C7) 



r^ )S (g) = l L{l qUB^ l^K + K 1 + K 2 )q 2 }, (CS) 



T s ^(q) = | A + ^ + q 2 + Ciq ^ ^ (CQ) 

r B)B (g) = ^L llQ 2 |s + (^K + K 1 + K 2 j g 2 } , (CIO) 

r tl t(g) = ^L {l q 2 jfi + (k + l -K^j g 2 } (Cll) 

r b)b (g) = ^-L {l q 2 (B + K iq 2 ), (C12) 
105 
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where A = \ + 1/0 and 5 = a /(p — a\. Although and 1^ contain a negative Co, both 
the growth rates increase monotonically with increasing q; therefore, a modulated pattern 
due to negative C will not appear. 

In the early stage, the concentration field and splay mode are coupled to each other via 
the off-diagonal terms. The eigenvalues of the coupled growth rates are T± = (T^ + r S)S ± 
y/iY^ — r SjS ) 2 + 4T Si ^T^ tS )/2. The peak positions indicated in Fig. [l] are calculated by 
solving dT x (q)/dq = 0, where T x (q) represents T + (g) ,r tt (g), and r bib (g). 

Shimada et al. derived similar dynamic equations for spinodal decomposition in an LCP 
solution [7] with the same initial equations; however, they derived their equations by ex- 
panding the Fokker-Planck equation directly and we obtain the free-energy functional and 
dynamic equations separately. However, it is noted that these differences are insignificant 
in the essential features of our results. 
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